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ABSTRACT. We prove two explicit bounds for the multiplicities of Steklov eigen- 
, values <Jk on compact surfaces with boundary. One of the bounds depends only 

■ on the genus of a surface and the index k of an eigenvalue, while the other de- 

pends also on the number of boundary components. Using spectral asymptotics 
for pseudo-differential operators, we also show that on any given smooth Rie- 
mannian surface with boundary, the multiplicities of Steklov eigenvalues cx^ are 
uniformly bounded in k. 



1. Introduction and main results 

1.1. Multiplicity bounds for Laplace eigenvalues. Let M be a smooth closed 
surface. For a Riemannian metric g on M we denote by 

o = Ao(s)<MsK---^teK--- 

the eigenvalues of the Laplace operator — A g . A classical result by Cheng in Q 
says that the multiplicities m^g) of these eigenvalues are bounded by quantities 
depending on the genus y of M only This result has been sharpened by Besson Q 
for orientable surfaces, and by Nadirashvili 11231 in the general case, to give the 
following version of the estimates for multiplicities: 

(1.1.1) m k (g) <2(2- x) + 2k+l, Jfc = 0,l,..., 

where % is the Euler-Poincare number of M. When M is a surface of zero genus 
these bounds are sharp for the first eigenvalue. 

The purpose of this paper is to prove analogues of such bounds for boundary 
value problems on surfaces. Such results seem to be missing in the literature, apart 
from the orientable zero genus case considered in ll23l 1X81 [T71 and Ofl, see also 
Remark fl.3.61 We are essentially concerned with the Steklov eigenvalue problem; 
problems with other boundary conditions are discussed at the end of the paper. 

1.2. Steklov eigenvalue problem. From now on let (M,g) be a compact Rie- 
mannian surface with a non-empty boundary. For a given bounded non-negative 
function p on the boundary dM the Steklov eigenvalue problem is stated as: 

(1.2.1) Au = inM, and ^- = cpu ondM, 
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where V is an outward normal. Denote by /I an absolutely continuous measure on 
the boundary dM with the density p, that is d\i = pds g . The real numbers a for 
which the solution above exists are eigenvalues of the Dirichlet form J \Vu\ 2 dVol g 
in the space L2(M,/i). Its spectrum is non-negative and discrete, see 0, and we 
denote by 

= a (g,n) < ax(g,n) ^...a k (g,n) 

the corresponding eigenvalues. This eigenvalue problem was considered in 1902 
by Steklov and since then was studied extensively; we refer to |2] and the recent 
papers Ifl4l l9l for a comprehensive list of references on the subject. 

I. 3. Main results. Our main result is the following theorem. 

Theorem 1.3.1. Let (M,g) be a compact Riemannian surface with a non-empty 
boundary. Let }X be an absolutely continuous measure on dM whose density is 
bounded. Then the multiplicity m k (g,}l) of the Steklov eigenvalue o k (g,}X) satisfies 
the inequalities 

(1.3.2) m k (g,il)^ 2(2-%) +2k + l, 

(1-3.3) m k {g,n)^2{2-x) + 2l + k, 

where % = % + l, and X and I stand for the Euler-Poincare number and the number 
of the boundary components ofM respectively. 

Note that % in the theorem depends on the genus y of M only. More precisely, 
it equals 2 — 2y or 2 — y with the respect to the cases M is orientable or not. Both 
inequalities (11.3.21) and (11.3.31 ) are similar in a way to the Besson-Nadirashvili mul- 
tiplicity bounds for Laplace eigenvalues on closed surfaces. The right hand-side 
of (11.3.21 ) is the same function of the genus of M as in (ll.l.ll ). This bound does 
not depend on any boundary data and, as we show, holds for other boundary value 
problems, see section [5] The second equality can be re- written in the form 

m k (g,n) ^ 2(2 -x) + k; 

it is specific to the Steklov problem, and for k > 21 — 1 is sharper than the first 
one. The proofs of both inequalities are built on the ideas due to l23l [18 ] and 
use the properties of nodal graphs. In comparison with other classical boundary 
value problems there is an additional difficulty related to the fact that for a Steklov 
eigenfunction there is no known local model for the nodal set at the boundary 
points. However, we show that, unlike for general harmonic functions, the nodal 
graph of a Steklov eigenfunction is still finite. 

For an annulus inequality (11.3.21) yields mi(g,n) 3; this bound is sharp, and 
the equality is attained on a "critical catenoid" constructed in In general, The- 
orem ll.3.T1 does not give sharp multiplicity bounds. It is an interesting question to 
understand whether m k (g,IJ.) is uniformly bounded in all parameters; see section 

II. 41 More specifically, one may ask the following question, cf. [15. Question 1.8]: 

Question 1.3.4. Does there exist a sequence of surfaces (M n ,g n ) with boundary 
measures ]X n such that mi(g n ,/J, n ) — > oo as n — > +oo ? 
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If such a sequence exists, by inequality (11.3.21 ) the corresponding genera y n of 
M n tend to infinity. Note also that the answer to an analogous question for the 
multiplicity of the first Laplace eigenvalue is positive Q. 

Remark 1.3.5. For a disk, Alessandrini and Magnanini proved in [1] the bound 
^ 2k, which is sharp for the first eigenvalue. In comparison, our result 
gives the bound mk(g,n) ^k + 2. Moreover, it can be slightly improved even 
further: for an even k we show that mk(g,n) ^ k+ 1, see Remark l3.3.2l 

Remark 1.3.6. While the present paper was at the final stage of preparation, a dif- 
ferent proof of inequality (11.3.21) for orientable surfaces appeared in lfT0ll20ll . The 
approaches behind all the proofs go back to the ideas of Cheng and Besson. At the 
same time, our proof that the nodal graph is finite is different from the one in ifTOTl : it 
is based on a topological argument and uses only Courant's nodal domain theorem. 
Besides, it applies to general boundary measures /x, see Lemma [3. 1.1 l and the dis- 
cussion in section [5] Note also that for non-orientable surfaces, inequality (11.3.21) 
is sharper than the estimate in lfT0ll20l . 

1.4. Asymptotic bounds for Steklov eigenvalues. Suppose that the boundary 
dM is smooth, and the weight function p in (11.2.11 ) is smooth and strictly positive. 
Then the Steklov eigenvalues can be viewed as the eigenvalues of a self-adjoint 
elliptic pseudo-differential operator of the first order; it sends a function on dM to 
the normal derivative of its harmonic extension multiplied by p ~ 1 . In particular, 
for p = 1, this pseudo-differential operator is precisely equal to the Dirichlet-to- 
Neumann operator on dQ.. Using Hormander's theorem on spectral asymptotics 
for pseudo-differential operators fl9l , we obtain the following result. 

Theorem 1.4.1. Let (M,g) be a compact Riemannian surface with a smooth bound- 
ary and /I be a measure on dM whose density p is smooth and strictly positive. 
Then the multiplicities mk(g,jl) of Steklov eigenvalues are uniformly bounded in k, 
i.e. there exists a constant C g n, depending on a metric g and a measure \l, such 
that 

m k (g,n) < C gjM for all A; = 0,1,.... 

The version of this result for Laplace eigenvalues is well-known, see ifTTl . When 
M is a disk, Theorem 1 1.4. II can be strengthened to the following statement. 

Proposition 1.4.2. Under the hypotheses of Theorem \1.4.1\ suppose that M is 
homeomorphic to a disk. Then there exists an integer K g n > 0, depending on a 
metric g and a measure jX, such that mfc(g, jtl) ^ 2 for all k > K g n. 

Note that the inequality above is sharp; it is attained on the unit Euclidean disk. 
The proof of Proposition 11.4.21 uses the unformisation theorem and the known 
asymptotics I24l l8l for the Steklov eigenvalues of the unit Euclidean disk. 

2. Preliminaries 

2.1 . Variational principle and Courant's nodal domain theorem. We start with 
recalling a variational setting for the Steklov eigenvalue problems. Given a Rie- 
mannian surface (M,g) and a measure /I on its boundary, the Steklov eigenvalues 
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can be defined by the min-max principle 



o*(&M) = inf SU P R s( M 'A0> 

A * +1 «eA* +1 



where the inrimum is taken over all (k + 1 ) -dimensional subspaces 1 C L 2 (M, /X ) 
formed by C°°-smooth functions, the supremum is over non-trivial u € A^ +1 , and 
Rg(«,ju) stands for the Rayleigh quotient 



Mention that the natural space of test-functions for the Rayleigh quotient above is 



here the second space in the intersection is formed by distributions whose deriva- 
tives are in Z,2(M, /i), see ll22l . The Steklov eigenfunctions can be then regarded as 
solutions of the equation 



understood as an integral identity, where the test-function <p ranges in Jzf. The 
equation above can be also viewed as a Schrodinger equation whose potential is a 
measure supported on the boundary of M. 

Let u be a Steklov eigenfunction. It is harmonic inside M, and, in particular, is 
C°°-smooth. By jV(u) we denote its nodal set, that is the set u~ l (0). Recall that a 
connected component of M\jV(u) is called its nodal domain. By maximum prin- 
ciple, it is straightforward to conclude that the closure of each nodal domain has 
a non-trivial intersection with the boundary dM. Further, by the strong maximum 
principle [12] any Steklov eigenfunction has different signs on adjacent nodal do- 
mains. For the sequel we need a version of Courant's nodal domain theorem for 
the Steklov eigenfunctions; see ||2T1 for the case of zero genus surfaces. 

Courant's nodal domain theorem. Let (M,g) be a smooth compact Riemannian 
surface with boundary, and /J. be an absolutely continuous Radon measure on dM 
whose density is bounded. Then each Steklov eigenfunction u corresponding to 
the eigenvalue <7k(g,n) has at most (k+l) nodal domains. 

Proof. Suppose the contrary. Then there exists more than (k+l) nodal domains 
G\ , . . . , Gk+i , Gk+2, ••• of an eigenfunction u. Consider the functions 



Since u is continuous and vanishes on dG\dM, it is straightforward to conclude 
that q>i € Jzf. Using it as a test-function in equation (12.1.2I ). we obtain that 




(2.1.1) 



Se = L 2 (M, pL ) n L\ (M, Volg) ; 



(2.1.2) 





(2.1.3) 




Now consider the sum 



k+l 

\jf = £ CCi<Pi. 



i=l 
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Since the closure of all nodal domains intersect dM, after an appropriate choice of 
the a,'s one can suppose that y is orthogonal in L%{M,ix) to all lower eigenfunc- 
tions. Then by definition of the eigenvalue Xk{pL,g) we obtain 

where in the last inequality we used relation (I2.1.3I ). Thus, we see that y has to 
be an eigenf unction itself. However, since it vanishes identically on Gt+2, it has 
to vanish identically everywhere by the unique continuation property fl2l . This 
contradiction finishes the proof. □ 

2.2. Local behaviour of harmonic functions; vanishing order. Let u be a har- 
monic function on M and x G M be an interior point. By the vanishing order of u 
at x we mean a non-negative integer, denoted by ord r (w), that is the order of the 
first non- vanishing derivative of u at x. The following statement is classical, see 
and lfl6l Theorem 4.1], and holds for solutions of rather general second order linear 
elliptic equations. 

Proposition 2.2.1. Let (M,g) be a smooth compact Riemannian surface with bound- 
ary, and u be a harmonic function on M. Then for any interior point xq G M there 
exists a non-trivial homogeneous harmonic polynomial P n of degree n = ord T0 (w), 
defined in a neighbourhood U ofxo, such that 

u(x) =P n (x-x Q ) + 0(\x-x Q \ n+1 ), 

where x G U. 

In the proposition above we assume that the neighbourhood U is such that the 
metric g\ v is conformally Euclidean. In particular, the property of being harmonic 
on U with respect to the metric g is equivalent to being harmonic with respect to 
the Euclidean metric. Now for a given positive integer I consider the set 

,jY\u) = {x G M | ord x (w) ^ £}. 

Using Prop. 12.2.11 in (51 Cheng shows that around a point xq G jV(u) the nodal 
set is diffeomorphic to the nodal set of the corresponding harmonic polynomial 
P n , which consists of n = ord Ao (u) lines meeting at the origin. In particular, the 
set jV 2 {u) consists of isolated points in the interior of M, and the complement 
jV {u)\jV («) is a collection of C°°-smooth arcs. Thus, the nodal set jV{u) can 
be viewed as a graph in the interior of M whose vertices are points x G jV 2 (u) and 
edges are connected components of jV x (u)\J/ 2 (u). In sequel we refer to Jfiu) 
as the nodal graph, meaning this graph structure. Mention that there are harmonic 
functions whose nodal graph is not finite. 

The following statement highlights the difference between Steklov eigenfunc- 
tions and arbitrary harmonic functions; it is a consequence of Lemma l3.1.1l in the 
next section. 

Proposition 2.2.2. Let (M,g) be a smooth compact Riemannian surface with bound- 
ary, and H be an absolutely continuous Radon measure on dM whose density is 
bounded. Then the nodal graph .jV{u) of a non-trivial Steklov eigenfunction u has 
a finite number of vertices and edges. 
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Mention that Prop. l2.2.2l does not bar possible pathologies of nodal set behaviour 
that do not occur for other classical boundary value problems. For example, it does 
not say that the closure of jV{u) intersects dM by a finite set. 

2.3. Graphs in surfaces: basic background. The purpose of this subsection is to 
introduce notation and collect a number of auxiliary facts used throughout the rest 
of the paper. Let S be a surface, possibly non-compact. Recall that a graph F C S 
is a collection of points, called vertices, and embedded intervals, called edges, 
such that the boundary of each edge belongs to the set of vertices. In addition, we 
assume that edges do not intersect and do not contain vertices. A graph is called 
compact if it is compact as a subset; it is called finite, if it has a finite number 
of vertices and edges. For example, for a non-trivial Steklov eigenfunction u the 
nodal graph jV{u), viewed as a subset in the interior of M, is not compact, since it 
contains edges approaching the boundary. 

Let r be a finite compact graph in S. For a vertex x G T its degree deg r (;e) is 
the number of edges incident to x; if there is an edge that starts and ends at x, then 
it counts twice. The number of edges e of a finite compact graph is given by the 
formula 

(2.3.1) 2e = £deg r (x), 

where the sum is taken over all vertices x € T. Connected components of S\F are 
called faces of T. The following inequality is a consequence of the standard Euler 
formula for a cell complex, see ifTTl p. 207]. 

The Euler inequality. Let T be a finite graph in a closed surface S, and v, e, and 
/ be the number of its vertices, edges, and faces respectively. Then the following 
inequality holds: 

(2.3.2) v - e + f>x , 

where % is the Euler-Poincare number of S. 

Finally, we recall the terminology for paths in graphs, which is used at the end of 
Sect. [3] By a path in a graph T we mean a continuous map : [0, 1] — > V such that 
0(0) and 0(1) are vertices, and if the image of intersects non-trivially with an 
edge, then it contains this edge. A path in T is called simple, if it has no repeated 
vertices and edges. A closed path in a finite graph is called the simple cycle, or 
circuit, if it has no repeated vertices and edges except for 0(0) = 0(1). A tree is 
a connected graph that has no circuits; its every two vertices can be joined by a 
simple path. 

3. Proof of Theorem I1.3.1I 

3.1. Reduced nodal graph. Let M be a smooth Riemannian surface with a non- 
empty boundary and M be a closed surface of the same genus, viewed as the image 
of M under collapsing its boundary components to points. By Jf{u) we denote 
the corresponding image of a nodal graph Jf{u); we call it the reduced nodal 
graph. More precisely, its edges are the same nodal arcs, and there are two types 
of vertices: vertices that correspond to the boundary components that contain limit 
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points of nodal lines, vertices of the first kind, and genuine vertices that correspond 
to the points in jY (k), vertices of the second kind. It is straightforward to see that 
the number of nodal domains of an eigenfunction u is precisely the number of the 
connected components of M\jY(u). Through out the paper we use the notation % 
for the Euler-Poincare number of M. It coincides with the quantity % + used in 
Theorem ll.3.11 and is called the reduced Euler-Poincare number of M. 

The following lemma is a basis for the proof of Theorem 11.3.11 It uses only 
Courant's nodal domain theorem, and holds for eigenfunctions of rather general 
boundary value problems. 

Lemma 3.1.1. Let (M,g) be a smooth compact Riemannian surface with bound- 
ary, and H be an absolutely continuous Radon measure on dM whose density is 
bounded. Then the reduced nodal graph JV (u) of a non-trivial Steklov eigenfunc- 
tion u is finite, i.e. it has a finite number of vertices and edges. 

Proof. Consider the reduced nodal graph ,jV(u) corresponding to a non-trivial 
Steklov eigenfunction u. To prove the lemma it is sufficient to rule out the oc- 
currence of: 

(i) infinite degree vertices of the first kind and 

(ii) the infinite number of vertices of the second kind 

in jV{u). We are going to construct new graphs in M by resolving vertices of 
JY{u) in the following fashion. Let x S <y¥ 2 {u) be a vertex of the second kind; its 
degree equals In, where n = oxA x {u). Let U be a coordinate ball centered at x that 
does not contain other vertices and such that nodal arcs incident to x intersect dll 
at 2n points precisely. We denote these points by y,-, where i = 0, . . . ,2n — 1, and 
assume that they are ordered consequently in the clockwise fashion. A new graph 
is obtained from Jfiu) by changing it inside U and removing possibly appeared 
edges without vertices. More precisely, we remove the nodal set inside U and 
round-off the edges on the boundary dU by non-intersecting arcs in U joining 
points y2j and yij+i- If there was an edge that starts and ends at x, then such a 
procedure may make it into a loop. If this occurs, then we remove this loop to 
obtain a genuine graph in M. It has one vertex less and at most as many faces as 
the original graph. 

Ruling out (i). Let us resolve each vertex of the second kind in j¥ (u) in the way 
described above. The result is a graph r in M whose only vertices are the first kind 
vertices in jV{u); we denote by v their number. Besides, it has at most as many 
faces as the reduced nodal graph, that is by Courant's nodal domain theorem at 
most k+ 1. Suppose that the reduced nodal graph has an infinite degree vertex of 
the first kind; then so does T. Let us remove all edges in T except for v + k + 2 — % 
of them to obtain a new finite graph, and denote by / the number of its faces. Since 
removing an edge does not increase the number of faces, we have / ^ k + 1. On 
the other hand, by the Euler inequality (12.3.21 ). we have 

f^e-v + x = k + 2. 

Thus, we arrive at a contradiction. 
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Ruling out (ii). Suppose the contrary; the situation described in (ii) occurs. Let v 
be a number of vertices of the first kind in ,jV(u). Let us resolve all vertices of the 
second kind except for v + k + 2 — % of them. The result is a finite graph V . Denote 
by V , e', and f the number of its vertices, edges, and faces respectively; then we 
have 

v' ^2v + k + 2-% and e ^ 2(v + k + 2 -%). 
Here in the second inequality we used formula (12.3. lb and the fact that the degree 
of each vertex x £ ,jV 2 {u) is at least 4. Combining the first two inequalities with 
the Euler inequality ( 12.3.21 ), we obtain 

/' ^e'-v' + x >k + 2. 
Thus, we arrive at a contradiction with the Courant's nodal domain theorem. □ 

3.2. Multiplicity bounds: the first inequality. We start with a lemma that gives 
a lower bound for the number of nodal domains via the vanishing order of points 
x € jV 2 {u). For the Dirichlet boundary problem on surfaces of zero genus it is 
proved in fl8l . 

Lemma 3.2.1. Let (M,g) be a smooth compact Riemannian surface with bound- 
ary, and /A be an absolutely continuous Radon measure on dM whose density is 
bounded. Then for any non-trivial Steklov eigenfunction u the number of its nodal 
domains is at least £(ord;t(w) — 1) +/£> where the sum is taken over all points in 
,jV 2 (u) and % is the reduced Euler-Poincare number ofM. 

Proof. Let ,J?(u) be a reduced nodal graph in M, and v, e, and / be the number of 
its vertices, edges, and faces respectively; by r we denote the number of first kind 
vertices. Using formula (12.3.11) . we get 

e ^ r + ^ord^w), 

where the sum is taken over x 6 jV 2 {u). Here we used the fact that the degree of 
each vertex of the first kind is at least two. Viewing v as the sum where the 

sum symbol is again over x € jV 2 (u), by the Euler inequality we obtain 

/ > e - v + £(ord x ( M ) - 1) + 
Since / is the precisely the number of nodal domains, we are done. □ 

The following lemma is a version of the statement due to [23 ]. 

Lemma 3.2.2. Let u\,... : U2 n be a collection of non-trivial linearly independent 
harmonic functions. Then for a given interior point x G M there exists a non-trivial 
linear combination £ a,-M,- whose vanishing order at the point x is at least n. 

Proof. Let V be the span of u\, . . . , ui n , and Vj be its subspace formed by functions 
u € V whose vanishing order at x is at least i, ord Y (w) ^ i. Clearly, the subspaces 
Vi form a nested sequence, V, + i C Vj. The statement of the lemma says that V n is 
non-trivial. Suppose the contrary, that is V n = {0}. Then the dimension of V is 
given by the formula 

dimV = 1 + £ dim(V-/V;- + i). 
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By Proposition 12.2. 1 1 the factor-space Vj/V;+i can be identified with a subspace of 
homogeneous harmonic polynomials of order i. In polar coordinates on R 2 such 
polynomials have the form 

P;(rcos 0,rsin 6) = ar'cos(nd) + br' sin(ra9), 

and in particular, form a space of dimension two. Thus, we obtain 

dimV <l+2(n-l)=2n-l. 

This is a contradiction with the hypotheses of the lemma. □ 

Now we prove the first inequality in Theorem ll.3.11 

m k {g,pL) < 2(2-2) + 2*+l. 

Suppose the contrary to its statement. Then there exists at least 2(2 — %) + 2k + 
2 linearly independent eigenf unctions corresponding to the eigenvalue Ou{g,pi). 
Pick an interior point i£M. By Lemma 13.2.21 there exists a new eigenfunction 
u whose vanishing order at the point x is at least 2 — % +k+ 1. Combining this 
with Lemma 13.2.11 we conclude that the number of the nodal domains of u is 
at least k + 2. Thus, we arrive at a contradiction with Courant's nodal domains 
theorem. □ 

3.3. Multiplicity bounds: the second inequality. The proof of the second in- 
equality is based on the following version of Lemma l3.2.1l 

Lemma 3.3.1. Let (M,g) be a smooth compact Riemannian surface with bound- 
ary, and H be an absolutely continuous Radon measure on dM whose density is 
bounded. Then for any non-trivial Steklov eigenfunction u the number of its nodal 
domains is at least 

max{2ord. v («) + 2% -21-2 \ x £ ^V 2 (u)}, 

where % is the reduced Euler-Poincare number of M and I is the number of the 
boundary components. 

We proceed with demonstrating the second inequality in Theorem [T3TTJ 

m k (g,n) ^2(2-x)+2l + k. 

Suppose the contrary to its statement. Then there exists at least 2(2 — x) + 21 +k+l 
linearly independent eigenf unctions corresponding to the eigenvalue Ou{g,pi). Pick 
an interior point x G M. By Lemma l3.2.2l fhere exists a new eigenfunction u whose 
vanishing order at the point x is at least 2 — # + /+[(& + l)/2], where the brackets 
stand for the integer part. Using the estimate in Lemma 13.3- 11 we see that the 
number of the nodal domains of u is at least k + 2. Thus, we arrive at a contradiction 
with Courant's nodal domains theorem. □ 

Remark 3.3.2. It is clear from the proof that actually a slightly improved inequality 
holds for the Steklov eigenvalue multiplicities: 

m k (g,n) <2(2- x) + 2l + k- [k} 2 , 



10 



MIKHAIL KARPUKHIN, GERASIM KOKAREV, AND IOSIF POLTEROVICH 



where \r1\2 is a function that takes values or 1 with respect to the cases n is odd 
or even respectively. 

The rest of the section is devoted to the proof of Lemma l3.3.1l We start with the 
following lemma, which is specific to the Steklov eigenvalue problem. 

Lemma 3.3.3. Let u be a non-trivial Steklov eigenfunction and x € ^V 2 (u) be a 
vertex in its nodal graph. Further, let T\ be a subgraph of jV(u) that is the union 
of all circuits in the connected component of x and all simple paths joining x and 
the vertices of these circuits. Then the degree ofx in T\ is at most 21 + 2 — 2%. 

Proof. Let vi, e\, and f\ be the number of vertices, edges, and faces of Ti respec- 
tively. Since every vertex in Y\, different from x, belongs either to a circuit or the 
interior of a simple path, its degree in Ti is at least 2. Thus, by formula (12.3.11 ) we 
have 

(3.3.4) 2 ei ^deg ri (x)+2(v!-l). 

Recall that every nodal domain of u has a non-trivial arc on the boundary. Since a 
face of Ti contains the union of nodal domains, it contains at least one boundary 
component of M. Since two faces can not contain the same boundary component, 
we have f\ ^ /. Viewing Y\ as a subgraph in the reduced nodal graph ^V(u), we 
can apply the Euler inequality to obtain 

e\ Kvi+fi-X^n + l-%- 

Now the statement follows by the combination of this inequality with relation (13.3.41) . 

□ 

Proof of Lemma \3.3. 1 1 Let x £ ^V 2 (u) be a vertex in the nodal graph, and F\ be a 
subgraph from Lemma 13.3.31 Further, let r 2 be a subgraph formed by all simple 
paths in the nodal set starting from x and approaching the boundary dM that do not 
intersect Ti except for x. This subgraph does not contain any circuits, and hence, 
is a tree. We denote by V2 and e2 the number of its vertices and edges respectively. 
Clearly, any nodal edge incident to x belongs either to F\ or T2, that is 

(3.3.5) 2ord^(w) = deg Fi (x) + deg rj (x) . 

We claim that the number of edges in F2 that are not incident to x is at least one 
less than the number of vertices, 

(3.3.6) e 2 - deg r , (x) > v 2 - 1. 

Indeed, this follows from the fact that T 2 is a tree, and that edges approaching the 
boundary have only one vertex. 

Finally, consider a subgraph T that is the union of Ti and r 2 . We use the notation 
v, e, and / for the number of its vertices, edges, and faces respectively. Clearly, we 
have 

e = ei+e2, v = Vi+V2 — 1. 
Combining these identities with relations (I3.3.4I )- (I3.3.6I ). we obtain 

e - v ^ 2ord. t («) - ]- deg Fi (x) - 1. 
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Using the bound for the degree from Lemma 13.3.31 we arrive at the relation 

e — v ^ 2ord A: (M) + % — I — 2. 

Finally, viewing T as a subgraph in the reduced nodal graph jV(u), we combine 
the last relation with the Euler inequality to obtain 

f^e-v + x-l>2ord x (u)+2x-2l-2. 

Since the number of faces / is not greater than the number of nodal domains, we 
are done. □ 

4. Asymptotic multiplicity bounds 

4. 1 . Proof of Theorem 1 1.4. 11 Let D be a "weighted" Dirichlet-to-Neumann oper- 
ator on dM defined as 

-\ * 

C°°{dM)3u^Du = p- 1 ^- £C°°(dM), 

dv 

where u denotes the unique harmonic extension of u into M. When p is smooth 
and positive, it defines a self-adjoint elliptic pseudo-differential operator of the 
first order whose eigenvalues coincide with the Steklov eigenvalues, see ll26l pp. 
37-38] and ll24l . Let N(X) be its eigenvalue counting function; it equals the num- 
ber of eigenvalues counted with multiplicity that is strictly less than a positive X . 
By the result of Hormander in |fl9ll25l . the function N(X) satisfies the following 
asymptotics (Weyl law): 

(4.1.1) tf(A) = A./ p(s)ds g + R(X), 

2% JdM 

where R(X) is a. bounded quantity in A > 0. Using this formula, we obtain 
m k {g,n) = YimN(X k + s)-N(X k ) = 



lim^- / p(s)dSg + R(X k + e)-R(X k )^2sup\R(X)\. 
e^o2n JdM 



Thus, the multiplicity m k {g : pL) is indeed bounded, and the theorem is demon- 
strated. □ 
It is interesting to know up to what extent the bound on m k (g,pL) depends on 
a metric and a boundary measure; in particular, whether there exists a universal 
constant (possibly depending on the genus of M) for which Theorem 1 1 .4. 1 1 holds. 

4.2. Proof of Proposition HA2l By the uniformisation theorem, we may assume 
that M is the unit disk and the metric g on M is conformal to the Euclidean met- 
ric gEuc- Since the Dirichlet energy is conformally invariant, from the variational 
principle it follows that the Steklov eigenvalies of (M,g) with a weight function p 
coincide with the Steklov eigenvalues of (M,gE uc ) with the a new weight function 
Po that depends on p and the values of g on dM only. By the results in 11241 181. the 
latter satisfy the following refinement of Weyl's asymptotic formula: 

2%k 2%k 

(4.2.1) Oat = 7 — n ru +°( k ~°°)i G2k+\ = - r — --rr+o(n. 
JdM Poinds J dM p Q {s)ds 
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as k — > oo. Thus, we conclude that for a large lc, the eigenvalues may be at most 
double. Note that for a disk all nonzero eigenvalues have multiplicity two, and 
therefore this result is sharp. □ 

Remark 4.2.2. The hypotheses of Proposition 1 1 .4. 21 on the smoothness of dM and 
p > are essential for the asymptotic formula (14.2.11 ) to hold. Even for domains 
with piecewise smooth boundaries the asymptotic properties of the spectrum may 
be quite different. In particular, by a direct computation one can show that for- 
mulas (14.2.11 ) fail for a square: for a large k the Steklov spectrum of a square is a 
union of quadruples of eigenvalues such that in each quadruple the eigenvalues are 
o(&~°°)-close lfT3l . However, no counterexample to Proposition 1 1 .4.2l is known for 
simply-connected surfaces with non-smooth boundaries, and it would be interest- 
ing to understand whether the result holds in this case as well. 

5. Other boundary value problems 

5.1. Eigenvalue problems with homogeneous boundary conditions. The method 
used to prove the first inequality in Theorem 1 1 . 3 . 1 1 uses only Courant's nodal do- 
main theorem and the behaviour of eigenfunctions in the interior of M; it largely 
disregards their behaviour on the boundary. The purpose of this section is to show 
that it applies to rather general boundary value problems. 

Let {M,g) be a compact Riemannian surface with boundary and L = (— A g ) + V 
be a Schrodinger operator, where V is a smooth potential. Denote by B a boundary 
differential operator of the form 

Bu = au + b—, 
av 

where a and b are smooth functions on dM and at least one of them does not have 
zeroes. We consider the following eigenvalue problem 

(5.1.1) Lu = Xu inM, and Bu = on dM; 

it is often referred to as the Robin boundary value problem, and in particular, con- 
tains the Dirichlet and Neumann problems as its special cases. By 

^ Ao ^ Ai < . . . fa < ■ ■ ■ 

we denote the corresponding eigenvalues. 

The following statement gives a bound for the eigenvalue multiplicities of prob- 
lem (15.1.11) that is independent of a Schrodinger operator L and, more interestingly, 
of a boundary operator B. 

Proposition 5.1.2. Let M be a smooth compact surface with a non-empty bound- 
ary. Then for any Schrodinger operator L and any Robin boundary operator B the 
multiplicity m^ of an eigenvalue A# corresponding to problem (15.1.11 ) satisfies the 
inequality 

(5.1.3) m k ^2(2-x) + 2k+l, 

where % = % + l, and X and I stand for the Euler-Poincare number and the number 
of the boundary components ofM respectively. 
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5.2. Details on the proof. We explain how the arguments and results in sections|2] 
and[3]carry over this general eigenvalue problem. First, Proposition 12. 2. 1 I holds for 
solutions of second order elliptic differential equations with smooth coefficients. 
In particular, it holds for eigenf unctions of problem (15- 1 - lb - Thus, the nodal set 
of an eigenfunction has a similar graph structure. These eigenfunctions also en- 
joy Courant's nodal domain theorem, see 0, and the arguments in section [3] show 
that their nodal graphs are finite. The version of Proposition 12.2.11 also shows that 
the statement of Lemma 13.2.21 holds for solutions of general second order ellip- 
tic equations, cf. |[23l Lemma 4]. The rest of the proof of the first inequality in 
Theorem 11.3.1 1 carries over without changes. 

Finally, mention that inequality (15 - 1 -3b holds also for eigenvalue problems with 
mixed boundary conditions. In addtition, one can also allow non-smooth bound- 
aries as long as the eigenvalue problem remains well-posed and Courant's nodal 
domain theorem holds. 
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